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IV. 

General Properties of certain Partial Differential Equations similar to those of 

Hydrodynamics. 

BY THOMAS CEAIG, Ph.D. 
U. S. Coast and Geodetic Survey, and Johns Hopkins University. 



Presented January 12, 1881. 



The following investigation was suggested by certain results which I obtained 
in a previous paper entitled " General Properties of the Equations of Steady Motion " 
which has been published by the United States Coast Survey. It was my intention 
in that paper to study steady motion only, but as I went on with the work, I per- 
ceived that many of the results were applicable to fluid motion in general, and that 
the whole subject admitted of a generalization which would present a number of 
interesting points. 

In Vol. LIV. of "Crelle's Journal" there is a paper by Clebsch entitled "fiber 
eine allgemeine Transformation der hydrodynamischen Gleichungen," from which I have 
taken the proofs of two properties of a certain determinant ; the proofs are also to be 
found in Baltzer's treatise and in Scott's treatise on determinants. I did not see 
Clebsch's paper until I had nearly completed mine, or I could have omitted one or 
two tedious processes that I have given by merely referring to Clebsch's article for 
the required demonstrations. It will perhaps be better, however, to leave this paper 
as it is, both because it is simpler and because the two papers are so unlike in their 
aims and methods of development, that, in order to merely refer to the article by 
Clebsch, it would be necessary for me to change my article throughout. Clebsch has 
in view merely a general transformation of the equations of hydrodynamics, while I 
have in view a general investigation of the " higher orders of motion " (vide Note at 
end) and of the properties of a certain set of Partial Differential Equations of the first 
order. 

The present paper is divided into two parts, of which Part I. is herewith com- 
municated to the Academy ; Part II. is not yet completed. In Part II. it is my 
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intention to generalize somewhat the methods and results of Part I. to make an 
application to the general theory of steady motion (n variables), and finally to apply 
all of the results of both parts to space of three dimensions, i. e. to the ordinary 
theory of hydrodynamics. 

(Just here I wish to call attention to the Note at the end of this paper; it 
properly belongs in Part II., but as I wish it made public at the earliest opportunity, 
1 have inserted it in Part I.) 

In Part I. an account is given of a set of partial differential equations similar to 
those of hydrodynamics, but containing n + 1 independent variables, or, if we may 
use the expression, the equations of fluid motion for space of n dimensions. 

Denote by x x , x 2 , x s , . . . . x n , t the independent variables of the problem ; by 
u u v^, . . . . u n , U, certain functions of these variables. For a problem in hydrody- 
namics the function U is given by 



u= V+ff, (1) 



V being the potential of external forces. The functions u are similar to the velocities 
in such an ordinary problem. Write 



.then 

Assume further 

and generally 

There will be of course 



2q = i h 2 + ii* + + Un 2 , (2) 

^ du x du 2 t du t du s 

&2 — dx a ~ ^j fi3 — 2xs ~ ck^' ete> > 

_duj _duz 
& dx k dxr ( 4 ) 



"j 



n (n — 1) 



of these quantities. The functions u are supposed to be connected by the relation 

du x d^ du n 

d& + ^ + ...- + ^ = 0. (5) 

Denoting now the operator 

— -4-— 4- ■ * 

by V' we nave > on observing that fa = — (■$, 
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~dx7Z ~ v M «- 



Between the three quantities f to , £ mn , int there exists for all values of I, m, n the 
relation 

WXj Wt^^j tttC^j \ / 

This relation holds throughout the three-dimensional space determined by the direc- 
tions x h x m , x n ; and as these are perfectly arbitrary, it is clear that a like relation 
holds in each of the Euclidean spaces determined by the directions x taken in groups 
of three. There does not appear to be any simple expression of this form connecting 
all of the ^-functions. 

Consider for a moment the differential expression 



Uydxx + Utflxz + . • . . + u„dx n ; 



(«) 



the conditions that this shall be an exact differential are (■■ 12 = £3 = = £„,„_i = : 

or, as these may all be written, 

d d d 

dx^ dx 2 ' ' ' ' ' dx n 



Mi, tl 2 , 



Un 



= 0. 



(8) 



If this quantity is an exact differential, say d^>, the function <f> must by virtue of (5) 
satisfy the differential equation 

V<£ = 0. 

If (a) can be made an exact differential by choosing a proper integrating factor, we 
must have 

Mlfji + U J&* + U *&J = ° ( 9 ) 

for all values of i, J, k ; or, as these may be written, 



1l\, M25 • • • . U n 

d d d 

dx^ dx 2 ' ' dx n 



u l ? u 2 j 



W» 



= 0. 



(10) 



SIMILAR TO THOSE OF HYDRODYNAMICS. 
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In a paper by Professor H. W. Lloyd Tanner "On the Transformation of a Linear 
Differential Expression," contained in the Quarterly Journal of Mathematics for 
December, 1878, it is proved that 



u x dx x + u 2 dx 2 + ....+ u n dx n 



is reducible to the form 



dvi + w 2 dv 2 + w z dv s + .... + Wjdv r 



if, and only if, 



OS) 



d 


d 


d 


dx^ 


cfej' 


' ' dx„ 


«i» 


U 2 , . . 


■ • u n 


d 


d 


d 


dx 1 ' 


dx 2 ' ' ' 


dx r 



= o, 



Ui, u 2 , . . . . u n 

where there are 2 r rows ; and that the same quantity is reducible to the form 

w i dv 1 + w 2 dv 2 + ....+ w r dv r , 



(11) 



if, and only if, 



u 



15 



u. 



2 J 



d_ d_ 



d 

dx n 



d_ d_ 



dx n 



= 0, 



(12) 



where there are 2 r + 1 rows. In both of these cases the functions w and v, though 
functions of x 1} x 2 , . . . . x n , are not mutually dependent. In the case in hand it will be 
shown that the expression (a) can always be reduced to the form (/3) when r = n + 1. 
The case when u x dx x + u 2 dx 2 + .... + u n dx^, is an exact differential d<f> gives for <f>, 
as already remarked, 

<P$ ,d% d 2 $ __ 

dx~? + dx~i + •••• + -fa} — 0. (13) 



I will give here, for convenience of reference, the transformation of this equation in 
polar coordinates. Some of the properties of the function $ will be found in another 

part of the paper. Denote by r, a u a^, <*„_! the polar coordinates of a point in 

space of n dimensions. Then write 

vol. x. 47 
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Xi = r cos a x> 

x 2 = r sin a x cos a 2 , 

x 3 = r sin a x sin a% cos 03, 



a^ = r sin a x sin 03 ... . sin a i _ 1 cos a* 



«„ = r sm oti sin 02 . . . . sin a n _ 1 . 
From these we have 

r = (^ 2 + X2 2 + . . . . + x n 2 f, 

X, 

COS Oh = - . 

1 r . 



cos a 2 



yV 2 - x? ' 
COS a 3 



03 ■ 

COS a 4 = 



sin 


<*2 
a 3 


= 


yV 2 - 

r 






sin 


Vr 2 - 


■ 2 - *i 2 


% 2 ) 


sin 


Vr 2 - 


W + 


*2 2 + %*) 




\/r 


! - (*i 2 


: + X,<) 


sin 


|/V - 


k = i 

- 3 
k = \ 

k — i — \ 


X k 
> 



4/V- S a* 2 i/V - S a* 2 



*=i 



tan (!„_! — 



x~ 



,5 



(14) 



V^-^ + a:/)' U3 ~ V H-fa-H-g,') ' ( 15 ) 



Now we have at once (writing for the moment a = r) 

2 . _ ^ v ^Vi ^t v feYj. _L o *± v *"* da * _L O ^ V *5» *i! _l 

V 9 — doj * \dx k ) "T" da^ ^ ^y T . . . . "I- ^ da(ldai 2, dXk dXk -f <* ^„7,;7 2, ^; ^ -f . . . 



1 datfla 2 dx k dx k 

+ Z v 2 «o + ~± v 2 «i + i: V 2 «2 + • • • • (16) 



Jao V «0 + &1 V «1 + d£ 



The summations in this formula are for all values of k from 1 to i. By performing 
the necessary operations we find quite readily 






dx k J a 2 sin 2 a x 



? \dx k ) ~~ <io 2 n sin 2 a,' J ~ *> 2 > * - 1 ' 



SIMILAR TO THOSE OF HYDRODYNAMICS. 367 

da, da i+1 . 

? &;^r~ ' * = o,i,2,....n-i, 

„ re — 1 2 _ (re — 2) cos a x 2 _ (re — 3) cos c^ 

V a ° ~~ ~^~> V a i — ^ sin - , V a 2 — ao l »sma,8in a a I ' 

<* _ l> - (i + 1)] cos a, - 2 _ 

V a J — ~ — T^i^i » V ««-i ~ U - 

sin a, II sin 2 a,- 
i = i 

Substituting in yf<f>, replacing cto by r, and multiplying through by r 2 we have 

^ 9 — r dr* "*" dB,.! -T S i n 2« rf 2 I" .... I" , = ;_! rf 2 I" .... I" , 2 ^ 

II sin 2 «,- ' II sin 2 «,- 

. ^/ _..<£<£ (re— 2) cos «j . d<j> (re — 3)cosa2 , , ^ [re— (t'+l)] cosa 4 

-+" T -j— In — 1) + -: : r* -7— — -. —„ [-.... + -; ,■ , . 

dr ^ ' da, sin a, aa, smitsiii'ii, <fe. • J= li~ • ■> 

1 ' * x ' sin«j- II snra,- 

+ . . . . '"* (16) 

or, arranging a little differently, 

r d*(r<f,) 1_ d_ / . d$ \ 1 ±_ ( • ^\ , 

|_ r <fr 2 "+" sin « x ^ Kl ^ sm °l do, J + sin « 2 . sin 2 «! da 2 \ Sm ^ da 2 J ~*~ ' ' ' ' 

sin a, n sin 2 «,. aa ' \ dcti J n sin 2 «,- da — 1 J 



+ [(n - 3) r ^ + (n - 3) 5^ ^ + (» - 4) 



COS « 2 1 C?<fr | 

J7 "1 



sin a 2 ' sm a i da 2 

i-in i*-r^;j sina .i-,-1 ^^••••^ S in« •'■?r*-« ■ d «»-J ^ ; 

* n Sin 2 «,• fein«„_ 2 jj gm 2 „ n ;_| 

Denote by T the region of all values of x t , x 2 . . . . jc„ for which 

^ 2 + x 2 2 + . . . . + x„ 2 = i? 2 (18) 

i. e. i? is the superior limit of r. Denote by S the boundary of T; then 8 signifies 
the totality of all systems of values for which 

Xi "T" 3?2 1" .... "" r" X n — -flj. I - '-"/ 

The equation y 2 ^> = holds throughout the entire region T and right up to the 
boundary 8. Then it can be proved just as in the case of three variables that there 
is one and only one real function <£ of the variables x, that for the whole region T is 
single valued and continuous, and at the boundary 8 coincides with an arbitrary 
function spread out continuously over 8. The first and second differential coefficients 
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of <f> will be continuous, and the latter will satisfy the equation y 2 <£ = 0. Denote by 
dT the element of the region T; then 

jfc==n— 2 n — k — 1 

dT — dx t dx 2 dx 3 .... dx n = r' i_1 ( II sin a da k ) d<t> n ~ x dr. (20) 

The element of the boundary 8 is dS, and this is readily found to be given by 

dS= R 1 - 1 Cn~ 2 sm~a rl da k ) da^^ (21) 

{R, a\, a 2 , . . • • a' n _! being the limits of r, ai, ag, . . . . a„_j) 
Introducing now an angle 8 defined by 

i = n — 1 j = k lc = n — l 

cos 8 = cos a a cos ai + 2 {cos a k + 1 cos c4 +1 II sin a 5 sin aj} + II sin a k sin a* (22) 
The introduction of this angle 8 enables us to obtain quite easily 

^ %/■■../- *<*->" , (23) 

2^ H ,?-! B (IP + r 12 -2 Br cos S)* 

<^! denoting the arbitrarily chosen surface value of <£. 

Further consideration of this function ^» will be deferred for the present, as a full 
investigation of its properties would require the development of the theory of 
spherical harmonics for a space of n dimensions. Equation (17) will afford the means 
of obtaining some of the more elementary of these properties. 

We will resume here the consideration of the quantities g jk . Observe that if we 
denote the quantities of the left-hand sides of equations (6) by a single letter, say by 
Vi> V*> • • ' ' ^ na ^ we have identically 

£ + £ + •■■•+£=«■ (24) 

If now from the i/s we form quantities similar to the £'s and call them {;$, £$ .... 
we shall have a set of equations of the same form as (6) with £ replaced by £ (1) and 
u replaced by rj. Calling the left-hand members of this new set rj£>, t$\ . . . ., and 
again performing the same operations, we get equations similar to (6) with the £ 
replaced by £ (2) , and for the right-hand members v/ 22 w, \7 22 u 2 , etc., where y 22 means 
that the operation v/ 2 has been performed twice. This process can of course go on 
indefinitely or until we come to such values of the ^-functions that 

ir^dxi + rfi > dx % +.... + rj^dXi 

is an exact differential. Instead of treating these functions u x , u 2 . . . . w 4 altogether 



SIMILAR TO THOSE OF HYDRODYNAMICS. 369 

we may consider them combined in groups of three each ; i. e. consider the motion in 

n (n — 1) (w — 2) 
6 

flat spaces of three dimensions instead of the motion in one space of n dimensions. 
Equation (7) shows the relation that exists between the rotation components g jk in 
such a space ; of course this relation holds for all of the Euclidean spaces which can 
be formed from the space of n dimensions, and therefore a similar relation exists for 
the n-dimensional space formed by suming the expression on the left-hand side of 
(7) for all values of I, m, n. Write for brevity 

n (n - 1) Q - 2) _ 

6 — 1 ' 

If the motion in any one of the / three dimensional spaces is rotational, it is of 
course rotational for the space of n dimensions. 
For the space containing x h x m , x n , we have 

lm _ $Tm» _ _±_ ( dn, ■ dn m . dlt H \ _ , 

Vi dx n dx m \dXi ' dx m dx n J V "i 

d (dn x dn m . dn n \ „ . K 

&T„ \ta, + s- + dx- n ) ~ V"n (25) 

dx, \dx, ^ dx m ~t~ dx n J V u i 

where v 2 is written for 

d 2 d* d* 

dx? * dxj ' dx/ 

For ordinary space the quantities in parentheses vanish and these reduce to well- 
known forms. * Write them in the forms 

dB „ 

Vi = dx, ~V u t 

Vr* = £ ~ V 2 «- ( 26 ) 

dB 



dilm 

dxi 


dx n 


dtmn 

ClX m 


__ d&a 
dx t 


djnl 

dx n 


_ d&n 

dx m 



y«-dx-~ v 2 «» 



Form now the quantities 



These operations cause the quantity E to disappear and leave U8 
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Vr»=V 2 £m (28) 

/ „2<? 

Perform the same operations on the if functions, and call the results rj" ', 77^, ■>?». 

// 2 dR 2 2 

< = v 2 £-v 2 - 2 ^ ( 29 ) 

// 2 dR „ , 

^ = v ^:~v ■«»• 

This operation can be carried on indefinitely, and we shall have in each case of an 
rj with an odd number of accents that the quantity R, or the expansion, as we may 
call it, disappears, while for the t/s with an even number of accents it appears in the 
form 

V dx 

k being here an odd number equal to the number of accents diminished by unity. 
For an incompressible fluid we have therefore a set of very simple formulas giving 
the relations between these various orders of motion. 

We shall now introduce the equations of "motion ; these are 

d U du x du t duy , . dic L 

dx\ ~ It + Ul dx~ t + " 2 dx~ 2 + *" Un dx~ n 

dTJ du 2 du 2 du 2 . du 2 

tot = * + lh d£ + *» ^ + ' • * + M « &„ ( 3 °) 

dTJ du n du„ du^ # , du^ 

dx~ n ~ Ht "*" Ul dx~ ± + Xh dx] + ' * ,_t " Un dx n ' 

Introducing the quantities £ and q and writing TJ — q = W, these equations may 

be written 

dW __du Y £ i . . . 4_ £ 



<&c 2 <fe 



-sz =-ji — %&2 + «8&aH t-««L ( 31 ) 



d TF du n t ^ _ .. 

^" — "^ WliXn M 2fc2» " * M„_ig: B _i >n . 
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du 



For steady motion the terms -^ vanish, and then multiplying these equations by 
u x , it 2 . . . u n respectively and adding, we have 



dW . dW , . dW 



dx x 



€&{£„ 



Equations (31) may be written in the form 

dW du k , 

= — - 4- H 

dx k dt ' ~ k 



(32) 



(33) 



Then, by taking equation (5) into account, we have for the determination of W 

(34) 






Take now any two of equations (10), say those indicated by the subscripts^' and k ; 
differentiate the /-equation for x k and the /^-equation for x j} and subtract the second 
result from the first; continue the process until all of equations (10) have been com- 



n(n— 1) 



bined two at a time by this process ; we thus arrive at a group of — t> 
of the form 



equations 



d$u . d£ u d$ n d£ n 

~dT + u ^'dx~ 1 + u ^dx~ i ^ r u n ^ 



+ 



t f ^ Ul _l_ ^ M A I t dug du t 

«2 \dx l "T" dX 2 J "T" W3 dx 2 "T" S14 <fo 2 "T 



+ £ 






523 TfoT ~~ fc24 fa. + 



(?«! 



rfXj 



£ — "=0 



or generally 



"3F + M i dx\ + W 2 ^ + 



d£jk _,_ _. <%& 
du x , u du<i 



+ W, 



d£ jk _|_ * /<&*/ , du t \ 
» <&„ "l" *& [dxj ' dx k ) 



t du x du 2 



du„ 



(35) 



A — "-0 



Since £# = £& = and ^ = — £%, we may write this equation in the form 



dt 






r=l 



Zrh 


£rj 


du r 

dx k 


du r 
dxj 



= 



(36) 



These equations give the vortex motion in the w-dimensional space ; for n = 3, 
they obviously become the ordinary equations for vortex motion. 
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The equations giving the motion of a fluid particle are, of course, 



~dt ~ u i>~dt =W2 - ■ -~dt 



= u„ 



(37) 



Writing for brevity 
equations (30) become 



D rf _i_ ^_i_ *^_i_ ^ 

dt. ~ df. + u l d^. ~r ^2 ^T + ' ' ' u n dx„ 



dt dt 



dU Bu x 



dxi 



dx 2 

d IT Du % 
dx„~ ~dt> etc * 



dx 1 ' dt 
Multiplying these by dx x , dx 2 . . . and adding, we have readily 

dU -\- dq — -j t (u^Xi + u 2 dx 2 + • • • u n dx„). 

Integrating from a point a to another yS, we have 

diJo. (ihdXi + u^dx^ + • • • u n dx n ) = \U + q]t 

If a coincide with fi, i. e. if we integrate around a closed curve, 

— y{ii x dx x + u 2 dx 2 + • • ' + u n dx n ) = 0, 



(38) 



(39) 



(40) 



(41) 



a theorem corresponding to that in hydrodynamics which states that the circulation in 
any closed curve moving with the fluid is independent of the time. As the principal 
object of this paper is to carry on and generalize the investigations begun in the 
paper on steady motion previously referred to, I shall leave for the present any 
further examination of the properties of the quantities dealt with so far and begin 
from a rather more general point of view. In steady motion a certain Jacobian ap- 
peared which entered into all of the work in a curious manner, this Jacobian was 
found from the functions 1} 2 , S , i. e. 



M x 


d0 1 


de x 


dx 


dy 


dz 


de 2 


d0 2 


d%. 


dx 


' dy ' 


dz 


de 3 


dOs 


de s 


dx 


' dy 


dz 



<M> 2 U/Vq U*l/ 2 j- 



(42) 



In the present more general problem we may introduce a set of functions 

7 di,0 a . .. . n -\ 



SIMILAR TO THOSE OF HYDRODYNAMICS. 



373 



of Xi , x 2 . . . . x n , t, denote the Jacobian of these by J, i. e. 



d6_ d0_ 
dO x d6y 
d6 n _ x d6 n _i 



' dx n 
dA 

(lX n 

dx„. 



= J. 



(43) 



The minors corresponding to elements in the first row may be denoted by 

-"4-n , -4i2 ,...., A ln , 
and in general for the minors to the elements in the i th row 

Of course we have now 

The minors -4 satisfy identically the equations 

dA,, dA ii . . dA,„ 

ll^ + -d^ + + ~dx^ ~ °> etc - 

If we substitute J. u ^ . ... A ln for ii 1} u 2 .... u n in equations (38) we 

at once 

1TT % n DA, j 
dU=Z ST dx J 



or 



j = n 



dA 



n j = n 



dA, 



dU= X~df dx J + 2 2 4i* -^f <&>,•. 



;=i *=i j=i 

Introducing the quantity 2 q which is now 

I= -«- ii + A-12 + .... + A yn 



(44) 

(45) 
obtain 

(46) 
(47) 



(47) becomes 



d(TT~q)=Z-df dxj+ 2 2^i 



i=i 



'j < *-l ^4 ---'-li 
* = 1 i=l 



d 

dx k 


d 

dXj 


Am 


A lj 



iXJbA 



For brevity write 





d d 


k = n 


— _— «— 


Kj = 2 ^it 


dx k dXj 


* = 1 


A-iic-Aij 



(48) 



(49) 



Assume a set of functions n — 1 in number 



vol. x. 



4>i , <E> 2 , 0„_, 

48 
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and write 

__ d$! dd 2 . , d6 n 

K i = *i &:+ ** ^-+ • . • • + $„_! ^, 

*i -*i^.+ ^d^ + .... + S>„_i ^, (50) 

Multiply each of these equations by A n , A n , . . . . A ln respectively and add: this 

gives 

A n K 1 + AuK 2 + .... + A ln K n = (51) 

a necessary condition, as is easily seen ; for if we multiply each K as given in (49) 
by the minor A having the same member for its second subscript and add, there 
results 

A n K x + A 12 K 2 + ....+ A ln K n =% £A u A lt 0gf - -£) . (52) 

An interchange of j with k in this should not alter the value of the right-hand 
member of the equation, but it obviously does alter it to the extent of changing its 
sign, and the quantity is therefore equal zero, or 

A 11 K 1 + A 12 K 2 + ....A ln K n =0. 

For the determination of the functions <E> we may proceed as follows. Denote the 
minors of A lh with respect to the element in the i + 1 th row and k th column by Af h ; 
then the quantity 

i = « — 1 , Q 

represents the determinant found from A lh by changing the column of elements con- 
taining differential coefficients with respect to % into a column containing differential 
coefficients with respect to x y In general then this determinant must vanish, as 
there will exist in it two columns containing the derivatives of the functions with 
respect to x y If j = h the determinant becomes simply A lh . In case j = h the 
determinant contains no derivatives with respect to Xj and for that case 

.A dx. A ™ = ~~ Alk 

since 

dA u __ _ (?i lt | 

dxk dxj 
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In all cases, then, when the above determinant does not vanish, it is equal to either 

dA 

A lh or — A lk . Multiplying it then by -^r and summing for h and k, from h = 1 to 
h = u and from k = 1 to k = « we have 

*Sl ft = l i^l <*"* dX J ~ ftfl dx i til " dx » ' 

Now the function Ej may be written in the form 

Kj = *S ^ u ' ^ ~~ a?i ^ A ^ ^ 

This differs from the second member of the last equation merely in sign. This equa- 
tion may be written in the form 

» = « — 1 (JQ, lk=n h=n 



t = l 

Now since 



dO- ( k — n h=n J A 1 



dxi ' dx 2 ' ' ' ' dx n 

the quantity in brackets can be written 

ZdV k )Z A lh AfA (y) 

k=i K h= 1 ) 

Now as we have for the determination of the functions <I> a set of equations of the 

form 

d9 l d0 2 dO n % n (dA xj dA lk \ 

^ ^ + ** &, + • • • • + *» ^ = fc fi ^ * (^ - ^7 J • 

Now observing the equations (a) and (/?) we have at once by comparing the coef- 

„.„ de- 
ficients of ^ 

^ = -2 2f^ (53) 



£i ifti dx * 
or, on taking (y) into account 

~:=n J rhz=n -\ 



k=n /f ( k = n 

<*>. = _ 

k 

Now since 



-"Lift 



<^1* 



1A .7 dSi 

a -r~ 

<lx k 



it is obvious that the quantity in brackets is the sum of the differential coefficients of 
2 q with respect to ^- and therefore finally 
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d dq d dq _i_ ^ dq 

~ ®i — & 1 ^ + ^ 2 ^ + + ^^7- (55) 

dx x dx 2 dx n 

If we multiply equations (50) by Ai_ ltl , -*4i_i, 2 , • • • •> -4»-i,» respectively and add, 
we shall obviously have 

*« = j{A i _ 1>1 K 1 + A i _^E 2 + ....+ Ai-^KJ (56) 

a form of <£* from which some interesting results may be obtained. Substituting^ 
for its value in (28) gives 

d(U-q) = 2^f dxj + 2 <hM> (57) 

and from this we obtain at once 

d(U-g) _ dAy dB 1 d6 2 de n _ x 

<£* — at +$ i(fe 1 + ^^ + — + *»-i "^T 

& 2 -^- + ^i^ + *2^ + .... + <*>„_! "3— (58) 



d(U-g) _dA ln i, , . ^, , , <»„■! 

tfo, ~ dt "•" ^i <fo n T" ^ &„ "+■ -r <P»-i ^ b 

For the case of steady motion the first terms on the right side of (58) drop out and 
the equations take the form 

d(U-g) _ , d9 x d0 2 <»„_, 

dx< =*i^ + ®2^ + .... + *„-i -^ (59) 

Multiplying by dx lf dx 2 , .... dx n respectively and adding, we see that for steady 
motion the functions $ are the differential coefficients with respect to 1} 2 , . . . ., 6 n _i 
of a function. 

w=u- q = w(6A.--.e n -i) 

i.e. 

(60) 





*, 


dW dW 

— dB x ' ^2 — M a • 


dW 


The operator 












^ n dx x ^ A ™ dx 2 ^ • ■ ■ 


+ A — 
• • ^ ^i» dx n 


is well known to be 


equivalent to 
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If we multiply now equations (58) by A n , A n , . . . ., A ln respectively we have by 
virtue of this relation between the two operators 

d(U-q) _dW _1 dq 

d6 >~ d$ >~ J dt \ bl ) 

and consequently 

W=U- g = fj§d6 + V(0 l e t ....O tt _ 1 t) (62) 

^ being a perfectly arbitrary functional symbol. This might be expressed differently ; 
since is the only variable function we must have 

4 dXl + ^ dx * + • • ' • + dV n dx » = dd 

&x dXi + di 2 dx * + • • • • + ix n dx « = ° 



from which follows 



_ A n d$ 
dr = dug. 



Substituting the value of d6 thus obtained in the above equation we are enabled to 
express the integral in terms of x-y, % 2 , . , . ., x n . In order however to obtain the 
value of the integral in (62) it will be necessary to express x x x 2 . . . . x n in terms of 
0, di, 2 , • • • .j n -i' To this end consider the Jacobian 

A = «(M--t-T> . (M) 

d (XyXz x n ) V J 

dW dW 
The minors of A corresponding to -fa , ~fa ■ ■ ■ ■ are the same as the minors of J" corre- 
sponding to -fa , -fa . . . . The minor corresponding to ^ is 

dW .... , dW A .„ , , dW A . ,„,, 

*f^i + ^^i + ..-- + a^^ (64) 

which is again a Jacobian. Denote the principal minors of A corresponding to ele- 
ments in the first row by A n , A 12 , . . . ., A ln and by A 41 , A,- 2 , . . . ., A in corresponding 
to elements in the I th row. Now as we know 

dA n rfA,- 2 dA m 

^+^+---'+ &„ = °5 ( 65 ) 
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the quantity (64) is = A ift , therefore (65) becomes 

J^^^^h - ( 66 ) 

From equations (58) we have, however, 

dW dA u dd,. , _ «W, , , , <»._! 

^ = ^T+*i^ + *^ + * * " * + *—i ^T 

Substitution of this in (66) gives a series of terms of which the first is 

This is = for all values of i other than i = 1. Similarly for all the terms until the 
i th ; the coefficient of <E> 4 is then =A lh . Now this substitution gives before reduction 

0=22^[^(^^ 1 + ^4+----+^- 1 -^- l +-^)J (67) 
on performing the above described reduction this becomes 

2 25 (^i**«) + 2 2 55 ( -ar ^) = ° (68) 

ft— i ft=i ^=i \ / 

Since the first and second minors of J", i. e. A lh and A'{ h , all satisfy identically the 
equation of continuity, (68) may be written in the form 



d<f>! 
•^ii dX! 


+ A ^ 


+ 


A d ** 


^ ^12 ^ 


+ 


^11 dx 1 


1 1 yl ^'- : 
1 A n dx, 


l + 



+^i»£ l +®x=o 



-H"" + 4.£ +e 2 = o 



+ 4* -^r +®=o 



(69) 



in which 



ft=« j=n , , . 

e. = 2 2 4? * ib? (70) 

ft=i ^=i 



Again, by aid of (53) we can replace (69) by 

a ^ + a - 1 A \- A ^ -4- ^ -X 

^n dx x ^ A ™ dx 2 ^ ^ A mdx n ^ dt ~ A i 



•^11 tfej ^ ^12 <fo 2 ^ 


_l_ A d ® 2 4_ ^* 2 


= *2 






= *»-! 



(71) 



SIMILAR TO THOSE OF HYDRODYNAMICS. 379 

in which 

^^ dA lH dA? 



dx 1 

~dt ~ w i' 




-di — u 2 , etc. 


dt — An* 




~dt = ""12 > e * 0. 




de k 

dt 


= 



The modifications undergone by (69) and (71) in the case of steady motion are 
obvious and need no remark. 
The equations 

have been replaced by 

dau flee, 

('73) 
If we have 



then 6 k = const, is an integral of these equations, and similarly for any other of the 
functions 6. If 

then 

d& k _ . d$ k d® k , d& k . d& k v _ 

lu~ A id^ + A zdx^~^''''~i~dl^~i~~dT~ x x~ ° 
and in this case 

<£> k = const. 

is also an integral of (73). If all of the functions Q 3 are independent of t, we shall 

have 

0x = const., 2 = const 6 n _ x = const. 

as integrals of (73). Suppose in equations (10) that the functions u were 

Uq Uj U% . . . . U 2n 

then it is known (see article of Professor Tanner's above referred to) that 

IIqCIXq i~ U\CLX\ 1 " I l U* 2n (XX 2n 

can always be reduced to the form 

clM + L 1 dM 1 + L 2 dM 2 + • f- L n dM n 

and 

dM dM 1 dM n 

w * = ^ + Z i ^ + + L n ~dx~ k (74) 

In this case we should also have 

fc=2» , j=2n fc=2» . , , . 

4^ _ 2)= 2i4 + l| 2 (i* - dxVj (V% - u * dx j) ( 75 ) 

fc=0 j=Q fc=0 \ J V 
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Now clearly 

du k duj i 



O/Xj CtXfc 



The double sum in (75) thus becomes 



dL r 


dL r 


dxj 


dx k 


dM r dM r 



r=2n 



r-0 



P=2re 
p-2n 



U, 



r=0 



dL r ^ n dL r 



p dx a 



dM r 



. dxo ax f> 



P=0 
p=2« 



dM r 



y u — r y 

■^i w p dx ? ^ dx„ 



dx n 



P=0 



P=0 



and since 

this is 



dXj 

dt 



Uj 



r=2n ( 
= 2 \ dM ' 


D d 

dt dt 


— dL r 


D d 

dt dt 


L* r Ld r 




M r M r 



Substituting this in (75) we have, after some reductions, 



,r rT dM r 4" r dM r -\ r 4" ( dL r „ r dM r ,\ 

* \U-q—*- 2 A *~ J = ,1 {n? ***--* dL,) (76) 



The equations 



du 



du„ 



~Tt + u odx-„ + w i 



du 
dx x 



+ 



dU 
dx 

dTJ du x dui du x 

dx L ~~It + n o dx t + u i dx~! + 



• du 

~t~ u 2n fa 

, ditt 

' ^2« dxT 



dTJ du. 2n . du 2n du. 

dxZ n ~ If "+" n o ~d^l + u i ~~' 



and 



dx 



du . 6?% </m 2 

«fo ' dx 1 ~*~ /*•»•- ' 



dx x 



+ 



+ «2» 



rf« 2 

(/av 



(/sc- 



are thus replaced by a set of the form 



dL r 



dL r 



~dt + u o id + «! 



<Of r 



rfilf r 



+ "0 AT + Mi 



dL r 
dx 1 

dM r 
dx 1 



+ 



+ 



<SE2 M 



</£ r _ dO 

+ W 2« *^ _ _ </£, 



(77) 



(78) 



(79) 



(78) remaining unchanged. In (79) r has, of course, all values from to 2n, and fl 
is an arbitrary function of the L% M's, and t. For U we have now 
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tt— dM ° 4- ^? r dM " \.i ^ ^ 4. X t dM " 4- O /-am 

P_l A." \ P — 1 } 

Equations (59), which correspond to the case of steady motion, give us, on multiplica- 
tion by Ai_ ltl , A -1,2 ? • • • • A-i,n an( i addition of the results, 

1 / , dW , A dW S X A dW \ /0,X 

* l =j{A i _ lA -^ + 4_ li2 ^ + • • • • + 4_ 1>B ^-J (81) 

expressions similar to (56) for the determination of the functions <J>. 
We have seen that if we have 

f = (82) 

then k is an integral of equations (73) or of 

dt ~ U ^ dt ~ U *' " ' ' dt ~~ W « V° rf J 

In order now to obtain an idea of the nature of the functions 0,0 X . . . . # M _i, 
consider a function of the n + 1 variables a^ , x 2 . . . . x n>t , which denote by/, i. e. 

This function/ which satisfies say 

/ = const. 

sustains the same relation to the space of n dimensions as a function of the three 
variables x, y, z, and also of t sustains to ordinary space, so for want of a better term 
I shall speak of it as a hypersurface. Now suppose the hypersurface / to be composed 
always of the same particles, it must obviously satisfy the differential equation 

f+«if 1 + ^f 2 + ----+^l[ = (84) 

(84) has as we know n independent integrals which we may denote by 

0, $i, 6 2 . . ■ • n -i 
The equations 

= const., X = const., etc. 

are, therefore, the equations of moving hypersurfaces which contain at all times the 
same fluid particles. The particles which are at any time on the intersection of any 
n — 1 of these hypersurfaces will remain on that intersection. As before we denote 
the Jacobian of the functions 0, 1} 2 , . . . . n _ 1 by J. So far as it has been con- 
vol. x. 49 
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venient from considerations of symmetry to write 0, 1} 2 . . . . 0„_ 1 as our functions, 
for the same reason it will now be more convenient to write 1} 2 , 3 . . . . n . In all 
that precedes, then, the function is identical with the new 1} and in general t _i 
(old notation) = % (new notation). The minors of J corresponding to 

dh Mi de 1 

dXi dx 2 dx n 

are A n , A 12 . . . . A ln respectively; and similarly for the other rows. 
If we solve the equations 

~dT + w i ^ + M 2 ^ + • • • • + u n ^ - 



+ «iSt*S+--.. + U.£ = (85) 



d6 9 , d0 2 . d$ a , . d0 2 

d$ n dO n dO d6 n 

~di + * te t + «* ^ ; + ■••• + «» & = 



we find 



f/Ml ( A ^ n ^ <ft ^ 21 ^ ^ dt ^ \ 

^ = "jf ^ + f ^-+.... + f A 2 | (86) 

^ ~ | rf* ^ ll » "•" <ft • A 2»-r -t- rf< ^l»„^ 



Substituting these values of ^« in 



dx 1 dx 2 ' ' dx n 



and we find by the rule for differentiating a determinant 

dj dJ . . dj __ 

dx[ + %h dx~ 2 + • • • • + M » dx~ n ~ dt 



dJ dJ . . dj dj 

^ dl + «» d£ + •••• + «. £T = ~ Tt ( 87 ) 



or J" satisfies the equation 

It would be interesting to determine the effect of the operator j upon the various 
minors of J, but it is sufficient for the present purpose to find the effect produced 
upon the first minors by this operator. The method of doing this is, in the case of 
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A u , to erase the first of equations (85) and solve the remaining n — 1 equations for 
u 2 ,u 3 , . . . ., u„. We have after doing this 



V% = - j |_^ + «x &;J 4S+ [* +«t ^J ^1+.... + |_a+«i eJ 



/)21 
A»-l, 



Differentiate these with respect to x 2 , x 3 , . . . . x n and add and we readily find 



(89) 



and in like manner 



eft ~~ Al dx, "•" ^ <fo 2 T" • ■ • • t" -4 ^ 

«fc ~~~ Au dx^ A Udx 2 " r "•" ^ <&. 

-P^is _ , ^s , A du s , , . &h 

dt — A * dx, ~*~ A ™ dx, "^ "i" ^ ln <fe„ (90) 

* ~ Au dx, ~T~ A ™ dx. 2 ^ "*" ^ lre <fo„ 

If in these equations we substitute for yl u the other principal minors A 2i , A 3i , . . . . A ni 

(i having all values from i = 1 to i = u in all of these) the equations will still hold. 

I do not write the proof of this out in full, as it is extremely simple and yet rather 

tedious. 

If we combine the equations 9 = const, in groups of n — 1 each, we shall of 

course obtain the n lines representing the intersections of n hypersurfaces represented 

by the equations 6 = const. Call these lines s 1} s 2 , . . . . s„; and denote by ds 1} 

ds 2 , . . . . ds n the n sides of an element of fluid occupying the space between the 

hypersurfaces 

6\ = a.\, 2 == aa • • • • v n = o-n 

and 

$1 = <Xi + dai, 2 = 02 + ^Ct2 . • . . n = On + cZot^. 

We must have now 

«i = #1 (#1 , # 2 • * . . x n , t) 

(91) 
<*! + <?ai = $i (oik + Xucfex + x 2 + X 12 ds 2 ....«„ + X ln <fe„ , t) 

the quantities X ll5 X^, . . . . X lB denoting the 'direction cosines' of ds t . Similar equa- 
tions of course hold for a 2 , a s , . . . . a n . By Taylor's theorem we now obtain 
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da i ~ j X « fa\ + X i2 di 2 + • • • • + X m ^J *n 
<ta>i— jAa 5- + A22 ^ -t- + Ag, ^- j ds 2 



(92) 



The values of the X's are known to be 



X =^ 



M2 



where 



2q 



, etc. 



(93) 
(94) 



2q = A{ l + A\ 2 + .... + A\ n . 
Similar expressions hold for X^ (; = 2 , 3 . . . . n and i = 1 , 2 . . . . n). The quantity 
2 g- can be expressed by a determinant viz. 

C22 C23 . . . . & n 
c 32 c 33 • • • • C Sn 



2g = 



in which 



C»2 C»3 . . . . C 



(95) 



C 22 



- f \d Xj ) > <* - f \ dXj ) > 



etc. 



dO, dOz 



d0 2 d6 4 



(96) 



° 23 = ^ fa, fa, c ™ ~ 2, fa fa etc 

j taking all values from j = 1 to j = n; and c te = c H . Substituting the values of 
Xu, X 12 , . . . ., X 1(l from equations (93) in the first of equations (92) and we find at 
once 

(97) 



2gv s?«, <$i 



similar expressions are to be found for c?s 2 , rfs 3 , . . . . ds n . Denote the projections of 

ds! on the axis by 

d$\d$\ <&£>; 

then from what precedes 

Mi 5 M2) • • • • • M« 



d$ v 



t^dx 1 ^ + Wgefc^ + .... + u n dx^ = -^ 



<#, rfft, 



d0. 



<*»/ <&,''*' *' dx n 



d$>, d$^ M. 

dx x ' dx 2 > ' ' ' "' cfo;„ 



(98) 
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with similar expressions resulting from the projection of ds 2 , ds 3 , .... ds n on the 
axes. Comparing these equations with (39) and writing 

U+q=W (99) 

gives 

j t (Mu + M12 + + Mi») jr = dW (100) 

and by (88) 

D dW 

j t (Mn + Mb + + Mi») = J~de 1 (101) 

and in like manner 



^ (M»l + Mn2 + .... + U n A nn ) =Jffl^ 



Assume a function P such that 



then since 



w=^r (102) 



we have 



An ■£- + A a ~ + .... + A ln ± = J ^ etc. (103) 

Mu + Ml2 + ....+ M„A» = ^ j^ + «Al j 

Ma! + u 2 A^ + + u n A 2n = J J ^- + i// 2 J (104) 

Mnl + M»2 +-... + »»A» = ^ I ffl n + ^» [ 

The functions ip contain 1} 2 , . . . . n but not t. Solution of (104) gives 
d0, UP I de, UP ) d6 n I dP i 

bearing in mind that 

II A m || = J— 1 (105) 

These equations for v^, u^ .... u n may be written in the form 



386 



GENERAL PROPERTIES OP CERTAIN PARTIAL DIFFERENTIAL EQUATIONS 



_ dP d0 1 d0 2 d0 n 

«i-a: + ^s+^&: + — + «k ^r 



dP 



n dx, ^ V2 dXi 



d6 1 



d$ s 



dx. 



dx<. 



dx 1 



d6. 



«» - xr + ^1 xf + fa ■&.. ; + + <£» xr 



rfcc, 



(106) 



M. 



+ & cbL + & xr + + *. xr 



<£»„ ' Vi dx H ' V 2 dx n ' V» rf^ 

and from these follows 

iixdx-L + u 2 dx 2 + ....+ u n dx n = dP + ^dO^ + tjf 2 d0 2 + .... + \fj n dd n (107) 

This result may be obtained in another manner which will bring into evidence some 
interesting properties of the functions with which we are dealing. The method to be 
employed I shall merely indicate, as it is similar to that employed in obtaining 
equations (61) and (62). The original equations of motion are 



now 



and consequently 



Dm, 
~df 



dU 
dx x > 



Du, 
~dt 



■>+ ~ dx* etC ' 



d _d0 1 d d$ 2 d 
dx ( dx, d$j ' dx { d$ 2 * 



+ 



dQ, _rf 
dxi d0 n 



Du 1 
~dT 

Du 2 d6 x dU 



d0, dU , d0 3 dU , 
dx x d6 1 dx x dd 2 



dt 



dt^ dU 
dx 2 d0 x dx 2 d0 2 



.+ 
.+ 



d0» dU 
dx x dd n 

dO n dU 
dx, d0 n 



(108) 
(109) 



(110) 



m n _d$ 1 du d$ 2 du de^du 

dt ~ dx n d0 x """ dx' n d6 2 "*~ ' * ' ' + dx n d0„ 

Multiply these by A u , A 12 , . . . ., A ln respectively; then by A 21 , A m . . . ., A 2n re- 
spectively, etc. ; add in each case and we shall obtain 



Du, Du., 



Du„ 



dt ' dt '••••» dt 



dfy d0 2 
dx x ' dx 2 ' 



d$ 2 
*' dx a 



dx x ' dx 2 ' ' * * ' ' dx n 



T dU 
= J df 1 ' etG - 



(111) 



SIMILAR TO THOSE OF HYDRODYNAMICS. 



387 



Now we have found that 



•P4, 

dt ~ ^ u dx x ~t~ ^dx.^ "^ ""»» <£* 






Using this equation, and also equation (103) and we are enabled to write at once 



DA n _ du x 
dt ~ J d0 1 

DA n _ du 2 
dt ~° dj x 



(112) 



I>An __ jd%K 

dt d d0 x 



Multiplying these by ih, u 2 , . . . . u n respectively and adding the sum to (111) we 
obtain 





U\ j t<2 j • • 


• . , u n 




d$ 2 d0 2 


cW t 


D 


i.tlA--l WtAsn 


* ' dx, 


dt 


• 






&» &. 


d6, 




dx x ' dx 2 ' 


' ' ' dx 



— J 



dW 
d0 1 



there are of course n 
written down. Since 



— 1 other equations of this form which can be immediately 



DJ 

dt 



= 



DW 
we have, if we write P = -jf integrate and solve the resulting equations for 

u x , U2, . . • • it n the same expressions as those marked (106). 

Note. — I desire here to call attention to the paper published by Professor Eowland in the American 
Journal of Mathematics, Vol. III. No. 3, in which a full account of the so-called higher orders of motion 
is contained. (I may also mention that I am informed by Professor Rowland that he gave a brief account 
of these motions at the meeting of the Scientific Association of the Johns Hopkins University on 
October 6th). 

The existence of these motions was discovered at about the same time independently by Professor 
Eowland and myself. He was led to consider them while treating of the equations of electromagnetism, 
and their connection with fluid motion. I was led to similar results by the study of certain possible cases 
of steady motion in viscous fluids, an account of which will be found in the Journal of the Franklin 
Institute for October, 1880, and also in the Philosophical Magazine for November, 1880. 
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I had devised no simple means of determining the lower orders of motion in terms of the higher 
orders, but had merely derived the higher orders from the lower by obvious processes of differentiation. 
In Professor Rowland's paper an elegant solution of the much more difficult problem is given. I may 
further say that up to the time of my becoming acquainted with Professor Rowland's research I had 
found no physical interpretation of these quantities so that the entire credit for the development of the 
theory of these motions is due solely to Professor Rowland. In the present paper the quantities corre- 
sponding to the higher orders of motion are denoted by the letter t\ with different accents and subscripts. 
My work on the subject of the higher orders of motion is to be found in the paper quoted at the beginning 
of this article, " General Properties of the Equations of Steady Motion." I believe that Professor 
Rowland has made no use of the integrating factors which I have there introduced. This part of the 
subject I hope to go into quite fully in Part II. of the present paper. 

Washington, D. C. 
January 2, 1881. 



